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ABSTRACT 

We stack 4.6 Ms of high spectral resolution XMM-Newton Reflection Grating Spectrometer 
spectra from galaxy clusters, groups of galaxies and elliptical galaxies. For those objects with 
a central temperature of less than 1 keV, we detect O VII for the first time, with a probability 
of false detection of 2.5 x 1CT 4 . The flux ratio of the O VII to Fe XVII lines is 1 /4 to 1 /8 of 
the emission expected for isobaric radiative cooling in the absence of heating. There is either 
a process preventing cooling below 0.5 keV, anomalous O/Fe abundance ratios, absorbing 
material around the coolest X-ray emitting gas or non-radiative cooling taking place. The 
mean N VII emission line is strong in the sub-keV sample. As the ratio of the hydrogenic N and 
O lines is largely independent of temperature, we measure a mean N/O ratio of 4.0 ± 0.6 Z Q . 
Although the continuum around the C VI lines is difficult to measure we can similarly estimate 
that the C/O ratio is 0.9 ± 0.3 Z . 
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1 INTRODUCTION 

The intracluster medium (ICM) is the few x 10 6 to ~ 10 8 K X-ray 
emitting plasma which fills the potential wells of objects ranging in 
scale from galaxy clusters down to massive elliptical galaxies. The 
X-ray surface brightness profile of these objects is often steeply 
peaked (e.g. Ste wart et al.|19 84). The implied mean radiative cool- 
ing times are less than 1 Gyr. In the absence of heating 10s to 100s 
of solar masses per year should be cooling below ~ 10 6 K to form 
a cooling flow (Fabia njl994) . 

The relative strength of X-ray emission lines is a good probe of 
the temperature distribution of the ICM. The major spectral lines in 
the 5 — 38A spectral range can be observed in bright objects using 
the Reflection Grating Spectrometer (RGS) instruments on XMM- 
Newton. 

The weak or missing Fe XVII emission lines, indicating mate- 
rial around ~ 0.5 keV, in RGS spectra showed that there is much 
less cool gas present in clusters than expected from cooling in the 
absence of heating (e.g. |Peterson et al^2 001 , 2003 ). However, these 
lines have been detected in several objects using deep observations 
( |Sanders et al.|2008||2"009") . There can be a wide range of x 10- 15 
in temperature, but less than 10 per cent of the emission expected 
for radiative cooling is seen. Active galactic nuclei (AGN) in the 
cores of clusters are believed to be the mechanism by which energy 
is supplied to the ICM to combat a large fraction of the radiative 
cooling (for reviews see Peterson & Fabian 2006, McNamara & 
|Nulsen|2007) . 

Fig. [Jj shows the emissivity of emission lines from certain 
strong X-ray emitting ions as a function of temperature. To cre- 
ate this plot we summed the emissivity of strong lines in the RGS 
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Figure 1. Emissivity of lines in the RGS waveband from various strongly 
X-ray emitting ions as a function of temperature at Solar metallicity. 



spectral range calculated using SPEX 2.00.11 (Kaas tra & Mewe] 
2000). The emissivity in FeXVII emission peaks at a temperature 
of around 0.5 keV. To observe gas at lower temperatures requires 
the detection of O VII lines. No significant detections of O VII have 
been reported in clusters, groups or elliptical galaxies, implying 
there is little material below 0.5 keV temperature. 
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Table 1. Details of the objects in the different temperature bins. Shown are 
the number of objects in the bin, the average redshift and the total RGS1 
plus RGS2 exposures. 



Temperature range (keV) 


Number 


<z> 


Exposure (ks) 


0-0.7 


9 


0.0063 


1240 


0.7-1 


9 


0.012 


1329 


1-2 


8 


0.027 


1847 


2+ 


34 


0.13 


4831 



Very deep XMM observations of bright high-metallicity cool 
objects may detect the presence of Ovil emission lines. The 
calorimeters planned for ASTRO-H and 1XO should also be able 
to detect these lines. In the absence of any such observations at this 
time we conducted a stacking analysis of the best observations in 
the XMM -Newton archive. 

In this paper we assume the relative Solar abundances of An- 
|ders & Grevesselp989) . 



2 SPECTRAL STACKING 

We took the processed RGS spectra from our analysis of the line 
broadening in 62 galaxy clusters, groups and elliptical galaxies 
(Sanders et al. 2010). These spectra were extracted from 90 per 
cent of the point spread function (PSF), a roughly 50 arcsec wide 
strip across the centre of the object, and included 90 per cent of the 
pulse height distribution. We removed from the sample the two ob- 
jects, Klemola 44 and RX J1347.5-1 14, which showed significant 
line broadening. 

We note that the sample is not complete and was originally 
chosen to contain objects with strong emission lines and a peaked 
surface brightness profile. The results are likely to be relevant for 
this sort of object, but perhaps not the average cluster, group or 
elliptical galaxy. 

For each object we created a fluxed spectrum with RGS- 
FLUXER, combining the two first order and two second order spec- 
tra. A fluxed spectrum uses a fixed wavelength binning and is 
divided by the effective area of the telescope to be in units of 
photon cm -2 s _1 A _I . We used a bin size of O.OlA when fluxing 
the spectra. We background-subtracted the spectra using a back- 
ground extracted from beyond 98 per cent of the RGS PSF (as in 
|Sanders et al.|2010| >. The RGSFLUXER tool is not normally used in 
the quantitative analysis of spectra because it does not take account 
of the redistribution into neighbouring channels. For the low resolu- 
tion spectra we examine here and because we only fit the emission 
lines with simple Gaussian functions, this is not a serious issue. 

The wavelength bins of the fluxed spectra were corrected for 
redshift (note that we did not change the flux, which is still the 
flux in the observed frame). We rebinned and added these spectra 
to make a stacked spectrum with 0.04A resolution. For each output 
bin we divided by the total number of input bins and spectra which 
were added for that bin, to construct an average spectrum. 

We constructed stacked spectra for the clusters in four differ- 
ent temperature bins using the average temperature measured from 
the RGS spectra (values are listed in Sand ers et al.|2010) . Table [T] 
shows details about the temperature bins. The average spectra for 
the objects in the temperature bins are shown in Fig [2] with So- 
lar metallicity model spectra (computed using APEC; Smit h et aT] 
|2001[ > at a representative temperature and spectral resolution. The 
panels show that as the temperature increases we see increasingly 
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Figure 2. The mean spectra of objects divided into four temperature bins 
using the RGS-measured temperature. In each panel we plot model spectra 
with similar spectra resolution and temperature at Solar metallicity. 

higher ionization states of iron at shorter wavelengths. There are 
no obvious lines in the data that are not accounted for in the model 
spectra. We see no evidence for any lines that may come from, for 
example, dark matter decay (e.g. Abazajian et al.|200l) . 

The spectra for the cooler objects are line-rich. To examine 
the lines in more detail we combine the two lower temperature 
bins. The mean spectrum for the objects with temperatures less than 
1 keV is shown in Fig.|3] On the plot we label the major lines and in 
the lower panel show model spectra at a range of temperatures. In 
this stacked spectrum we detect the recombination and forbidden 
lines of the O VII triplet at 21.6 and 22. lA, respectively. O VII is a 
good indicator of material below 0.5 keV (Fig. [T|. In addition we 
find the Cvi lines at 33. 7 A, although this part of the spectrum is 
noisy and the continuum uncertain. C VI has been seen previously 
in M87 ( [Werner et al.|2006| >. 

Due to the rebinning and fluxing process we cannot propa- 
gate the error bars robustly from the original spectra. Instead, we 
used the standard deviation of the data points in the spectral region 
around around the O VII lines (20.8 to 24. 5A, excluding the lines 
themselves) to estimate the error bars on each spectral point. A con- 
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Figure 3. (Top panel) Mean of spectra of all objects below 1 keV temperature. This spectrum has a total XMM-Newton exposure time of 1.2Ms. (Bottom 
panel) Model spectra with a range of temperature computed with APEC at Solar metallicity. 



stant model fit over this region gives a y} of 108.4 with 92 degrees 
of freedom. If we add to the model two Gaussian components of 
the same width at the wavelengths of the Ovil lines, the x 2 re- 
duces to 87.4 with 89 degrees of freedom (a reduced x 2 of 0.98). 
An F-test suggests that the probability of of this improvement in 
fit happening by chance is 2.5 x 10~ 4 . The separate probability for 
the 2 1 .6A line is 1 .4 x 10~ 3 and 0.06 for the 22. 1 A line. Similarly, 
the probability of false detection of the C VI line is 0.02. 



3 DISCUSSION 

The detection of O VII emission shows that there is material below 
~ 0.5 keV in the X-ray waveband, at least in objects with a mean 
central temperature below 1 keV. The signal does not seem to be 
stronger when the sample is split into two bins (Fig.[2j, which im- 
plies that it is unlikely that there is a single dominating object or 
that only the cooler objects show O VII emission. 

The flux of the Ovil recombination and forbidden lines for 
the sub- 1 keV temperature sample, as obtained from the Gaussian 
model fit, are shown in Table [2] The FeXVII lines at 16.78 and 
17.08A were also fitted and the flux calculated (the second line is 
actually two unresolved lines at 17.055 and 17.1 A). For this fit, 
we fixed the continuum at the value obtained from the O VII lines 
due to the proximity of neighbouring lines. We did not measure 
the fluxes for the Fe XVII lines around 15A because it is difficult to 
separate the continuum and these lines. 



Table 2. Line fluxes observed in the mean sub-1 keV temperature spec- 
trum compared to cooling flow models (units are 10~ 5 photon cm~ 2 s -1 ). 
Model fluxes are for a 1 M Q yr cooling flow, cooling from 1.5 keV at So- 
lar metallicity at a redshift of 0.0 1 , absorbed with a Galactic column density 
of 10 20 cm~ 2 . 



Line 


Observed 


Model (APEC) 


Model (SPEX) 


Fe XVII 16.78 


2.52±0.07 


1.1 


1.8 


Fe XVII 17.08 


5.39±0.08 


2.5 


4.1 


O VII 21.6 


0.31 ±0.09 


1.2 


1.0 


O VII 22.1 


0.20 ±0.08 


0.63 


0.67 



We calculated a model spectrum where gas is cooling isobari- 
cally at a rate of 1 Mr.-, yr~ 1 at a redshift of 0.0 1 , from a temperature 
of 1.5 keV to 0.08 keV (assuming H = 70kms~ 1 Mpc" 1 ). The 
spectrum is absorbed photoelectrically at an equivalent Hydrogen 
column of 10 20 cm~ 2 , for both the APEC and SPEX models. The 
flux values for the same lines are shown in Table [2] For the APEC 
model the Fe XVII to O VII flux ratio is 2.0 (summing the pairs 
of lines) and 3.5 for the SPEX model. In our average spectrum the 
ratio is 16. Therefore if the FeXVII and Ovil emission is caused 
by cooling, assuming Solar metallicity ratios and no extra absorp- 
tion on the coolest gas, there is 1/4 to 1/8 of the emission below 
0.5 keV temperature than expected from the emission above those 
temperatures. 
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Table 3. Hydrogenic line fluxes observed in the mean 0—1 keV spec- 
trum compared to thermal models (units are 10~ 5 photon cm~ 2 s ). Model 
fluxes are for a Solar metallicity plasma at 0.5 keV, absorbed by 10 20 cm~ 2 , 
with an Xspec normalization of 10~ 3 and at a redshift of 0. 



Line 


Observed 


Model (APEC) 


Model (SPEX) 


O VIII 19.0 


2.62 ±0.08 


17 


17 


Nvn 24.8 


0.75±0.11 


1.2 


1.4 


CVI 33.7 


0.40 ±0.14 


2.3 


2.1 



There may be more material cooling to a low temperatures if 
the cooler gas is surrounded by absorbing material. The edge for 
oxygen absorption occurs over a wavelength range 22.6 to 23.7A 
(seen in the TBNEW absorption model of |Wilms et al.|2000") . This 
edge is at longer wavelengths than the O VII triplet. If the oxygen 
metallicity is low in the cooling regions, this may lead to weak 
Ovil emission. Oxygen lines are the main coolant in the RGS 
waveband but iron emission is substantial at longer wavelengths 
in the far UV. The abundances in the central parts of these objects 
is likely to be unusual, particularly because of depletion onto dust 
grains (e.g. Canning et al. 2010). Non-radiative cooling of the ICM 
may also lead to there being less emission from cooler tempera- 
tures than expected, for example, because of the mixing of X-ray 
emitting material with colder gas (Fa bian et al.|2 002). 

The average metallicities can also be estimated for the sam- 
ple. Note from Fig.[T|that the ratios between the hydrogenic O VIII, 
N VII and C VI emissivities are the same above a temperature of 
0.25 keV. If there is little material below this temperature (as is 
indicated by the lack of O VII emission) then the ratios of these 
lines should be good indicators of metallicity. The Mg XII hydro- 
genic lines around 8.4A could be used for abundance measurements 
above 1 keV temperature but they are unfortunately blended with 
FeXXIV. 

Table [3] lists the line strengths for the hydrogenic lines in the 
mean spectrum. The error bars on the spectral data points for the 
N VII and C VI lines were estimated from the standard deviation 
of the data points around each line. When measuring the flux in 
the the O VIII line we assumed the same continuum level that we 
measured for the O VII lines. We fitted separate continuum levels 
for N VII and C VI (which are lower by 3 and 20 per cent from the 

VII continuum level, respectively). 

We also tabulate in Table [3] the line strengths for a simulated 
0.5 keV spectrum at Solar metallicity, absorbed by a column den- 
sity of 10 2l) cm~ 2 . The APEC and SPEX models show similar fluxes 
for the lines. At Solar metallicity there should be a O/N flux ratio 
of ~ 14. We observe a ratio of 3.5, which implies that the average 
nitrogen abundance relative to oxygen is 4.O±O.6Z ; 0. We note that 
the O VIII emission line is resonantly scattered which could lead to 
an underestimation of the O/N ratio, although this effect is rela- 
tively small for an elliptical galaxy like NGC 4636 (Wern er" et al.| 
|2009l >. 

The N VII to O VIII line ratio appears to be largest in the 0.7- 

1 keV temperature range (Fig. [2](. The strength of the N VII line 
relative to O VIII decreases with increasing temperature. This could 
either be due to real abundance changes with temperature, or it may 
be due to the fact that we are sampling larger regions for objects at 
higher redshifts. 

We also applied spectral fitting to the combined spectrum for 
objects with temperatures of less than 1 keV. Although there are 
difficulties in interpreting the results from such a spectral fit, we 
confirmed that nitrogen to oxygen ratio is greater than 3 times the 



solar value. As the continuum at long wavelengths is uncertain, it 
is unclear whether O/H metallicities are low or whether the N/H 
values are high. The N/Fe ratio is around 1.7, although the FeXVII 
lines are poorly fit. The two temperature component model gives 
temperatures of 0.34 and 0.71 keV, with the emission measure of 
the lower temperature component 14 per cent of the emission mea- 
sure of the hotter component. 

Nitrogen abundances appear to vary considerably in detailed 
observations of individual objects, although the exact values are 
dependent on modelling. |Xu et al.|p002| l found that NGC 4636, a 
relatively cool elliptical galaxy at ~ 0.6 keV, does not have strong 
N VII emission, with N/Fe at approximately Solar ratios. |Werner| 
|et al.|{2009) , however, found 2.5 Z s N/Fe in the core of this ob- 
ject and high nitrogen abundances in four other elliptical galaxies. 
In Centaurus the N/Fe ratio is around 3Z© ( |Sanders et al.| 2008). 
Centaurus has a mean RGS temperature of 1.6 keV. The N/Fe ra- 
tio in Sersic 159 is 0.1 Z (de Plaa et al.||2"00"6| >, with a mean 3 
keV temperature. M87, which has a mean temperature of 1.8 keV 
shows a N/Fe ratio of 1.6Z (Werner et al^|2CjO"6) . The N/Fe ra- 
tios for Abell262, Abell3581 and HCG62 (relatively cool objects 
with RGS temperatures less than 1.5 keV), are close to Solar values 
(Sanders et al. 2010). It is unclear why there is such a wide range 
of N/Fe values. A future study of a sample of bright objects using 
the ratio of the N VII to O VIII lines would be less model dependent 
than using N/Fe. 

Possible other observational evidence for high nitrogen abun- 
dances are the high [N II] A6584/HG! ratios in the nebulae surround- 
ing some central cluster galaxies, which can be as much as 3 (Hatch 
|et al.|2007) . Ratios as high as this are difficult to achieve in models 
without super-Solar nitrogen metallicities. 

Examining the carbon metallicity, the C/O flux ratio implies 
an carbon to oxygen abundance ratio of 0.85 ± 0.32 Zq . 

Nitrogen and carbon metallicities are useful to measure be- 
cause these elements are not created by supernova explosions (un- 
like the other elements detected in the X-ray waveband) and are 
typically created in low to intermediate mass stars (for a discussion 
see Molla et al. 2006). Nitrogen can be made in a primary process in 
intermediate mass stars, where oxygen and carbon generated before 
CNO cycling is converted to nitrogen. This should give a nitrogen 
enrichment which is independent of the initial metallicity of the 
star. In addition, nitrogen can be created by a secondary process, 
where carbon and oxygen present at the birth of the star is trans- 
formed into nitrogen. This should give a nitrogen enrichment pro- 
portional to the initial abundance. Oxygen enrichment is thought to 
be mostly due to Type II supernovae, produced from massive stars. 

The high average value of N/O observed in the sub-keV tem- 
perature sample is at the extreme end of the distribution of H II 
regions and stars in late-type galaxies shown in |Henry et al.| ( [2000] >. 
However, N/Fe has been found to be slightly supersolar in a sample 
of early-type galaxies (Toloba et al. 2009). The high N/O values we 
observe may either be due to enhanced nitrogen abundances or low 
oxygen abundances. 



4 CONCLUSIONS 

We examine the mean RGS spectrum of objects from the XMM- 
Newton archive in different temperature bins. For those objects 
with a RGS-measured temperature of less than 1 keV we detect for 
the first time O VII emission. This is indicative of material below 
0.5 keV. 

If the gas in the objects is cooling without any heating, there 
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is 1/4 to 1/8 of the expected O VII emission compared to the 16.8 
and 17.1 A FeXVII flux. This may imply that a small fraction of 
the gas is cooling in these objects, the abundance ratios of O and 
Fe are anomalous, or that there is absorbing material around the 
coolest gas. An alternative explanation is non-radiative cooling of 
the coolest gas. 

The ratio of the hydrogenic lines in the spectrum are inde- 
pendent of temperature, allowing estimates of the mean N/O and 
C/O ratios. In the sample of objects below 1 keV the N/O ratio is 
4.0 ± 0.6 Z B . The C/O ratio is 0.9 ± 0.3 Z B . 
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